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Abstract. In turbulent dynamos the production of large-scale magnetic fields is accompanied by 
a separation of magnetic helicity in scale. The large- and small-scale parts increase in magnitude. 
The small-scale part can eventually work against the dynamo and quench it, especially at high 
magnetic Reynolds numbers. A one-dimensional mean-field model of a dynamo is presented 
where diffusive magnetic helicity fluxes within the domain are important. It turns out that this 
effect helps to alleviate the quenching. Here we show that internal magnetic helicity fluxes, even 
within one hemisphere, can be important for alleviating catastrophic quenching. 
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The magnetic fields of astrophysical bodies like stars and galaxies show strengths 
p_ _,, which are close to equipartition. The scale of the magnetic field is larger then the dis- 
sipation length and reaches the order of the size of the object. The mechanism which 
X/^ \ creates those fields is believed to be a dynamo. Large- and small-scale magnctic hclic- 
ity with opposite signs are created. For high magnetic Reynolds numbers, ReM, this 
makes the dynamo saturate only on a resistive time scale and reduces th e saturation 



field strength much below equipartition (jBrandenburg fc Subramanianll2005l ). This effect 



is called catastrophic quenching and increases with increasing Reynolds number, because 
for the Sun ReM = 10 9 and for gala xies ReM = 10 14 . This suggests that helicit y has 
to be shed. Observations have shown ( Manoharan et al. 1996 : Canfield et al.| [T999) that 



helical structures on the Sun's surface are more likely to erupt into coronal mass ejections 
[ (CMEs). This suggests that the Sun sheds magn etic helicity by itself. 



¡> ■ In our earlier work ( Brandenburg et all 120091 ) we have considered a one-dimensional 
mean-field model in the z-direction of a dynamo with wind-driven magnetic helicity flux 
■ where the wind increases with distance from the midplane. Magnetic helicity evolution 



is taken into account by using what is known as the "dynamical quenching" formalism 
that is described in our earlier paper and in refcrcnccs thcrcin. Wc augmcnt thcsc studics 
by imposing a constant shear throughout the domain which facilitates the growth of the 
magnetic energy. We perform simulations in one hemisphere where we set the magnetic 
\ field in the z-direction to be symmctric (S) or antisymmetric (A) at the midplane. The 
L* ■ outer boundaries are set to either vertical field (VF) or perfect conductor (PC). The free 
. [ parameters are the dynamo numbers C a and Cs- By varying both numbers we find the 
critical values for dynamo action (Fig.Q]and Fig.[5]). The critical values for C a decrease 
when the shear increases. This is expected, since larger shear leads to stronger toroidal 
- - ■ field which enhances the dynamo effect. 

In the rest of this paper we study in more detail the case Cs = —10 and consider 
positive values of C a . In Fig.[3] we compare the dynamical a quenching model (using 
a magnetic Reynolds number of ReM = 10 5 ) with the standard (non-catastrophic) a 
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quenching where a cx 1/(1 + B / B eq ) with B being the mean field and B eq the equipar- 
tition value. Note that in the former case, the energies cross. Nevertheless, the A solution 
is stable in both cases - at least for C a ^ 10. This is demonstrated in Fig.[4l where we 
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Figure 1. Critical values for the strength of 
the forcing C a and the shear Cs for which 
dynamo action occurs for the cases of verti- 
cal field boundary conditions and antisymmet- 
ric (solid, red) and symmetric (dashed, blue) 
equator. The circles and squares represent os- 
cillating and stationary solutions respectively. 
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Figure 2. Critical values for the strength of 
the forcing C a and the shear Cs for which 
dynamo action occurs for the cases of per- 
fect conductor boundary conditions and anti- 
symmetric (solid, red) and symmetric (dashed, 
blue) equator. The circles and squares rep- 
resent oscillating and stationary solutions re- 
spectively. 
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Figure 3. Comparison of the bifurcation diagrams for dynamical and standard a quenching. 
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Figure 4. Evolution of the magnetic energy (left) and parity (right) for a solution that was 
initially even about the midplane (S or quadrupolar solution), but this solution is unstable and 
developed an odd parity (A or dipolar solution). 



show that after about 40 diffusive times, = 40, where rj t is the turbulent magnetic 

diffusivity and k\ the basic wavenumber, the magnetic energy E decreases and the parity 
P swaps from +1 to —1; see iBrandenburg et al. for details on similar studies 
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Figure 5. Profiles of qk, «m, and their sum for the A and S solutions at C a = 10. 
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Figure 6. Time-averaged magnetic helicity fluxes, Fi and F m , of fluctuatingjind mean fields, 
for the A and S solutions, respectively, at C a = 10. Note that F¡ + F m ~ 0. 

The crossing of the energies in the dynamical quenching model is somewhat surprising. 
In order to understand this behavior, we need to look at the profiles of the a effect; see 
Fig.[SJ In this model, a is composed of a kinetic part, aj¿, and a magnetic part, %, 
which has typically the opposite sign, which leads to a reduction of a — a^ + om- The 
quenching can be alleviated by reducing om, for example when the divergence of the 
magnetic helicity flux of the small-scale field, F¡, becomes important. 

Naively, we would have expected that the A solution should have a larger energy, 
because only this solution allows a magnetic helicity flux through the equator; see Fig.|6l 
This is however not the case, which may have several reasons. Even though thc magnetic 
helicity flux flux small at the equator (z = 0), there can be significant contributions from 
within cach hcmisphere which contributes to alleviating the catastrophic quenching. Thc 
details of this will be address in more detail elsewhere. 



References 

Brandenburg, A., & Subramanian, K. 2005, AN, 326, 400 
Brandenburg, A., Krause, F., Meinel, R., et al. 1989, A&A, 213, 411 
Brandenburg, A., Candelaresi, S. & Chatterjee, P. 2009, MNRAS 398, 1414 
Manoharan, P. K., van Driel-Gesztelyi, L., Pick, M. & Demoulin, P. 1996, ApJ 468, L73 
Canfield, R. C, Hudson, H. S. & McKenzie, D. E. 1999, Geophys. Res. Lett. 26, 627 



